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ABSTRACT 

We derive the black hole mass density as a function of redshift from the hard X-ray AGN luminosity function 
assuming that massive black holes grow via accreting the circumnuclear gases. The derived black hole mass 
density matches the measured local black hole mass density at z = 0, which requires the average radiative 
efficiency of AGNs to be ~ 0.1-0.17. The massive black holes in most faint active galactic nuclei (AGNs) and 
even normal galaxies are still accreting gases, though their accretion rates are very low. Radiatively inefficient 
accretion flows (RIAFs) are supposed in these faint sources, which should radiate mostly in the hard X-ray 
band. We calculate the contribution to the X-ray background from both the bright AGNs and the RIAFs in 
faint AGNs/normal galaxies. Our calculations show that both the observed intensity and spectral shape of the 
XRB with an energy peak at ~ 30 keV can be well reproduced without including the emission of Compton- 
thick AGNs, if the massive black holes in faint AGNs/normal galaxies are spinning rapidly with a ~ 0.9 and 
accreting at rates m~1.0-3.0x 10~ 4 . It indicates that less than ~5 per cent of local massive black hole mass 
density was accreted during radiatively inefficient accretion phases, which is obviously only an upper limit, 
because Compton-thick AGNs have not been considered. If the same number of the Compton-thick AGNs 
with log Nh = 24-25 as those with log Nh = 23-24 is considered, the fraction of local black hole mass density 
accumulated during inefficient accretion phases should be lower than ~2 per cent. The constraints of the XRB 
can provide upper limits on the average accretion rate for inactive galaxies. 

Subject headings: galaxies: active — quasars: general — accretion, accretion disks — black hole physics; X-rays: 
diffuse background 



1. INTRODUCTION 

There is evidence that most nearby galaxies contain mas- 
sive black holes at their centers, and a tight correlation was 
revealed between central massive black hole mass and the 
velocity dispersion o f the galaxy (Ferrarese & Merrittl 120001 : 
iGebhardt et alJ uOOO). which strongly suggests co-evolution 
of massive black holes and their host galaxies. This correla- 
tion, together with the correlation between black hole mass 
and the host galactic bulge luminosity, are widely used to es- 
timate the black hole masses in galaxies. 

The growth of massi ve black hole s may probably be linked 
to accretion processes (Soltani fT982h . lYu & Tremaind (2002) 
estimated the black hole masses from the stellar velocity dis- 
persions of galaxies measured by the Sloan Digital Sky Sur- 
vey (SDSS) using the empirical relation between black hole 
mass and the velocity dispersion, and the local black hole 
mass density was derived. They further calculated the black 
hole mass density accreted during optical bright quasar phases 
using an optical quasar luminosity function (LF), and found 
that the accreted mass density is consistent with the local 
black hole mass density estimated from the velocity disper- 
sions, if a radiative efficiency ^0.1 is adopted for quasars. 
This implies that the growth of massive black holes may be 
dominantly through accretion during optically bright quasar 
phases. The optical quasar LF is directly linked to the ac- 
cretion history of massive black holes, however, it has over- 
looked many faint AGNs and type II AGNs. The X-ray lumi- 
nosity functions (XLFs) derived from the surveys in the soft 
X-ray band (< 3 keV) may have missed many obscured (type 
II) AGNs, while the hard X-ray surveys (~ 2-10 keV) can 
trace the whole AGN population includin g obscured type I I 
AGNs and extend to low X-ray luminosity (lUeda et ai1l2003l) . 
thus providing a useful tool to explore the accretion history of 
AGNs in the universe. Many works on XLFs were carried out 



in either the soft X-ray band (< 3 keV) (e . g..lMaccacaro et all 
[19911 iBovle et alJ[l99l iPage et al.lfl997t iMivaii et alJl2000h 
or the hard X-ray band (> 2 keV) (e.g.. lBoyle et alJll998h 
ICowie et al.ll2003llUeda et alj|2003l) . 

The cosmological X-ray background (XRB) is mostly con- 
tributed by AGNs (lHasingerlll99"8l: ISchmidt et al.l[l998h . In 
the most popular synthesis models of the XRB based on the 
unification schemes for AGNs, the cosmological XRB from 
< 2 keV to more than several hundred keV can be fairly well 
reproduced by using a series of template spectra of AGNs 
with different obscuration, which are derived by assuming 
an intrinsic AGN spectrum consisting of a power-law X- 
ray spectrum with a n exponential cutoff around several hun- 
dred keV (see, e.gjMatt & Fabianlll994t[Madau et alj[l994t 
I Comastri et all 119951. iGillietalJ 1 1 999fc lUeda et all l2003l) . 
Di Matteo & Fabian ( 1997) have alternatively proposed that 
the hard XRB above 10 keV may be dominated by the ther- 
mal bremsstrahlung emission from the advection dominated 
accretion flows (ADAFs) in low-lumin osity AGNs. Further 
detailed ADAF spectral calculations (Di M atteo et alJfl999h 
have shown that many sources at redshift z ~ 2 - 3 with 
ADAFs accreting at the rates close to the critical value would 
be required to reproduce the observed XRB spectral shape 
with an ene rgy peak at ~ 30 k eV. 

Recently. lUeda et all d2003l) derived a hard X-ray luminos- 
ity function (HXLF) from a highly complete AGN sample (2- 
10 keV), which includes both type I and type II AGNs (except 
Compton-thick AGNs). Based on this HXLF, their synthesis 
models can explain most of the observed XRB from the soft 
X-ray band to the hard X-ray band around several hundred 
keV. Their calculations slightly («10-20%) underestimate the 
relative shape of the XRB spectrum around its peak intensity 
(see Fig. 18 in their work). Such a discrepancy can be ex- 
plained provided the same number of Compton-thick AGNs 
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with logJV H = 24-25 as those with logA^ H = 23-24 is in- 
cluded. As the number density of Compton-thick AGNs is 
still quite uncertain, the residual XRB could be attributed to 
both the Compton-thick AGNs with logNu = 24-25 and the 
RIAFs in faint AGNs/normal galaxies. Although their rela- 
tive importance is still unclear, it is obvious that the spectral 
shape of the observed XRB cannot be solely attributed to the 
RIAFs/ ADAFs in Seyferts. 

There are a variety of studies exploring the evolution 
of AGNs base d on either optical qua s ar LFs or XLFs, 
or both (e.g.. lHaehnelt & Reel 119931; lHaiman & Menoul 
2000t Kauffma nn & Haehneltl l2000t lYu & Tremaind I2002F 
Wvithe & Loebl 120031: iMarconi et ail 12004 iHopkins et ail 
2005). A common conclusion is that the timescale of AGN 
activities is short compared with the Hubble timescale, though 
the quantitative results on the bright quasar lifetime vary from 
~ 10 7 to ~ 10 9 years for different investigations. The standard 
optically thick accretion disks are present in bright AGNs, 
as their accretion rates are high. The AGN activity may be 
swit ched off while t he gases near the black hole are exhausted 
(see|Naravan 20021 for a review, and the references therein). 
When the accretion rate m (m = M/Mndd) declines below a 
critica l value m cr i t ~ 0.01, the standard disk transits to a RIAF 
(e.g jNarayan & YiF l995). RIAFs are optically thin, very hot, 
and their sp ectra are pe aked in the hard X-ray band. In a pre- 
vious work. iCaol d2005l) calculated the contribution to the hard 
XRB from the RIAFs in the faint AGNs. It was found that the 
accretion rate m must decrease rapidly to far below the criti- 
cal rate m cr i t after the accretion mode transition, otherwise, the 
RIAFs may produce too much hard X-ray emission to match 
the observed hard XRB. This implies that the growth of mas- 
sive black holes d uring radiative l y ineffi cient accretion phases 
is not important. Hopkin set al.l (|2006) considered the distri- 
bution of local supermassive black hole Eddington ratios and 
accretion rates, and they found that black hole mass growth, 
both of the integrated mass density and the masses of most 
individual objects, are dominated by an earlier, radiatively ef- 
ficient, high accretion rate stage, and not by the radiatively 
inefficient low accretion rate phases. 

In this paper, we will explore quantitatively on how much 
mass in the local massive black holes was accreted during 
their radiatively inefficient accretion phases from the con- 
straint of the XRB. The cosmological parameters SIm = 0.3, 
ft\ = 0.7, and Hq = 70 km s" 1 Mpc" 1 have been adopted in 
this work. 

2. BLACK HOLE MASS DENSITY 

The HXLF given bv lUeda et al1J200l is so far most suit- 
able for our present investigation, as it includes both type I 
and type II AGNs (exce pt Compton - thick AGNs). The HXLF 
in 2-10 keV derived bv lUeda et all (120031) is 
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All the parameters for our present adopted cosmology are as 
follows: A = 5.04(±0.33) x 10" 6 Mpc" 3 , logL»(ergs s" 1 ) = 
43.94^1, 7l =0.86 ±0.15, 72 = 2.23 ±0.13, pi = 4.23 ± 
0.39; p2 = = -1.5, z c = 1.9, and logL„(ergs s " 1 ) = 44.6 
dUedaetal.1120031 All the sources described by lUeda et ail 
|2003|)'s HXLF have hard X-ray luminosities L 2 " 10 keV be- 
tween 10 41 5 - 10 48 ergs s" 1 . In this work, we refer to 
the sources described by the lUeda et ail d2003 )'s HXLF as 
active galaxies, while all the remainder with L^" 10 keV < 
10 415 ergs s" 1 as inactive galaxies. 

The cosmological evolution of black hole mass density 
caused by accretion during active galaxy phases is described 
by 



dp b A h fe) 
dz 
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where <£>(Lx,z) is the HXLF given by Eq. (UJ), p^ h (z) (in units 
of Mq Mpc" 3 ) is the black hole mass density accreted during 
active galaxy phases between z and Zmax, and e is the radiative 
efficiency for AGNs. The ratio / cor of the bolometric lumi- 
nosity to the X-ray luminosity in 2- 10 keV is given by 

log/ cor = 1.54 + 0.24£ + 0.012£ 2 -0.0015£ 3 , (6) 

where C = logLb Q i/Lp - 12 , and L^ Q \ is the bolometric lumi- 
nosity dMarconi et al.ll2004l) . 

The central massive black holes in inactive galaxies may 
still be accreting gases, though the accretion rates may be 
very low as their hard X-ray luminosities L^f 10 keV are be- 
low 10 4L5 ergs s" 1 . The Eddington accretion rate Mead = 
(M bh /M )L Q , Edd /?7 eff c 2 , where L Q , Edd = 1.38 x 10 38 ergs s~\ 
and a conventional radiative efficiency r/ e ff = 0.1 is adopted 
here. This radiative efficiency ?7 e ff is only for the definition of 
the Eddington rate ME dd , and it can be different from the real 
radiative efficiency e in Eq. (0. When the mass accretion rate 
m (rh = M/Medd) is below a cr itical value ffi c rit, the standard 
disk converts to a RIAF (e.g.. iNarayan & Yil[l995h . RIAFs 
should be present in most, if not all, of those inactive galaxies. 
Besides the black hole mass density p£ h (z) accreted in active 
galaxies, we need to calculate the black hole density p^ h (z) 
accreted in inactive galaxies. Assuming an average dimen- 
sionless mass accretion rate mf^ t for those inactive galaxies, 
we can calculate the black hole mass density p^ h (z) accreted 
during inactive galaxy phases between z and z max by 



bh(z) 



aver 
inact 



^O.Edd(l- 



,RIAF 



)df 



dz M 77 e ffC 2 dz' ^ 

where e RIAF is the radiative efficiency of the RIAFs in inac- 
tive galaxies, provided the black hole mass density of inactive 
galaxies p™ h act (z) (in units of Mq Mpc" 3 ) as a function of red- 
shift z is known. Here, p™£ ct (z) is the black hole mass den- 
sity of inactive galaxies at z, while p^ h (z) describes the mass 
density accreted during inactive galaxy phases between z and 
Zmax- The radiative efficiency g RIAF is usually much lower than 
that for standard thin disks. Our numerical calculations for 
RIAF spectra show that e RIAF is significantly lower than 0.01, 
so we approximate e RIAF ~ in Eq. (0 while calculating the 
black hole mass density. 

The total black hole mass density p^{z) accumulated 
through accretion both in active and inactive galaxies is: 



P^(z) = p£ h (z) + pZ h (z), 



(8) 
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i.e., a sum of black hole mass densities accreted in active 
and inactive galaxy phases. In order to calculate p^(z), we 
need to know the black hole mass density of inactive galax- 
ies jO{, n h act (z) as a function of z (see Eq. |7j- Assuming the 
growth of massive bla ck holes is dominated by accretion (e.g., 
lYu & Tremaind l2002h . the total black hole mass density is 
given by 

Pbh(z) ^ Pbh(z) + Pbh(Zmax) = Pbh(Z) + Pbbfe) + Pbhfemax), (9) 

where pbh(Zmax) is the total black hole mass density at Zmax- 
The black hole mass density for active galaxies in the co- 
moving space at redshift z can be calculated from the XLF 
by 



fraction of local black hole mass density was accumulated in 
inactive galaxies (~16 per cent, for m™Zt = 0-001, see Table 
1). If mg^ t < 3 x 10" 4 , less than ~ 5 per cent of local black 
hole mass density was accreted in inactive galaxies. At low 
redshifts with z < 2, the black hole mass density is dominated 
by the black holes in inactive galaxies (see Fig. |2). In Fig. [2] 
we find that the derived black hole mass densities p™ ct (z) in 
inactive galaxies at low redshifts mainly depend on the value 
of the measured local black hole mass density p^ 31 adopted. 

5r 
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/corlx ^^ dLogLx M Mpc- 
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dLogL x 



(10) 

where <&(L x ,z) is the HXLF given by Eq. O, m™f is the av- 
erage dimensionless accretion rate for the active galaxies. The 
black hole mass density for inactive galaxies jO b n h act (z) can be 
derived by subtracting that for active galaxies p^(z) from the 
total black hole density pah(z). The black hole mass density 
of inactive galaxies /O b n h act (z) can therefore be calculated by 

PbTk) = Pbh(z) - (&Z) = pUz) + f&fc) + Pbh(Zmax) " P&Z). 

(ID 

Assuming all gal axies are active at z = Zmax (e-g-, 
iMarconi et al.ll200 4). the black hole mass density pbh(Zmax) = 
Pbhfemax) is available by using Eq. ([Toll. At z = Zmax, p£ h = 
Pbh = 0- Thus, we integrate Eqs. (0 and Q simultaneously 
over z from z max using Eqs. (TToT > and ( fTTT i. and the black hole 
mass densities as functions of z are finally available provided 
the three parameters: m a ™ r , tnf^ t , an( l e are specified. We 
adopt Zmax = 5 and m a ™ r = 0.1 in all our calculations. The 
value of the radiative efficiency e is tuned to let the derived 
total black hole mass density pbh(z) match p b °^ al at z = 0, for a 
given m a ™ t . 

The local black hole mass density p^ 3 ' was estimated from 
the stellar velocity dispersion measured with the SDSS by us- 
ing the relation between black hole mass and the stellar veloc- 
ity dispersion. The resulting local black hole mass densities 
vary from ~ 2 x 10 s to 5 x 1 5 M(n Mpc" 3 (e.g.JSalucci et al 



19991 1 Yu & Tremaine 2002; S heth et al.! 2003; Marco ni et al 



2004). In some previous works, the black hole mass densities 
as functions of redshift were derived to match the measured 
local black hole mass density by integratin g AGN LFs (e.g., 
lYu & Tremainel l2002t IMarconi et all |2004|) . The difference 
between these previous works and ours is that we have in- 
cluded the contribution to the local black hole mass density 
from accretion in inactive galaxies. 

Our calculations are carried out for two different val- 
ues of the local black hole mass density: p b ° cal = 2.5 x 
10 5 M Q Mpc -3 , and 4.0 x 10 5 Mq Mpc" 3 , respectively, con- 
sidering the different values of of? 2 ' derived in the literature 



(e.g., ISalucci et all 11999: 


Yu & Tremainel 120021 ISheth et all 


2003; Marconi etal J 12004 


). We find that the radiative effi- 



ciencies e ~ 0.07-0.12 are required, so that the total black 
hole mass densities accumulated through accretion can match 
the measured local black hole mass densities (see Figs. Q] 
and Table 1 for a summary of the results). The derived radia- 
tive efficiency depends only weakly on m™!£ v because most 
of black hole mass was accreted during active galaxy phases. 
If the average mass accretion rate is high, a significant 




FIG. 1 . — The total black hole mass densities pbh(z) accumulated through 
accretion as functions of redshift z (solid lines), while the black hole mass 
densities accreted during active galaxy phases only [p^ h (z) + Pbhfenax)] are 
plotted as dashed lines. The upper set of curves is derived for piS^ = 
4.OxlO 5 M Mpc" 3 , while the lower one is for p b ° h cil1 = 2.5 X 10 5 M Q Mpc" 3 . 
The different colors correspond to different models: black (Model A), blue 
(Model B), red (Model C), green (Model D), and yellow (Model E). The de- 
scription of different models is summarized in Table 1 . 




FIG. 2. — The black hole mass densities 



p b n h act (z) for inactive galax- 
ies as functions of redshift z (solid lines), while the dashed lines represent 
the black hole mass densities p bb (z) in active galaxies. The upper set of 
curves is derived for pjjjf 3 ' = 4.0 X 10 5 Mq Mpc -3 , while the lower one is for 
^bh 31 = ^ x i " 5 Mq Mpc -3 . The different colors correspond to different 
models: black (Model A), blue (Model B), red (Model C), green (Model D), 
and yellow (Model E). The description of different models is summarized in 
Table 1. 



3. RIAF SPECTRA 

The RIAF spectrum is almost proportional to black hole 
mass Mbh for a given dimensionless mass accretion rate m. In 
order to calculate the contribution of the RIAFs in all inac- 
tive galaxies to the XRB, we need to use the X-ray spectrum 
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lx(E) of a RIAF around a 10 M© black hole accreting at the 
rate »i?™ t as a template spectr um. We employ the approach 
suggested by Manmoto (2000) to calculate the global struc- 
ture of an accretion flow surrounding a massive black hole 
in the general relativistic frame, which allows us to calculate 
the structure of an accretion flow surrounding either a spin- 
ning or non-spinning black hole. All the radiation processes 
are included i n the calculation s of the global accretion flow 
structure (see Manmoto 2000, for details and the references 
therein). However, the va lues of some par ameters adopted 
are differ ent from th ose in iManmotol d2000) (see discussion 
in §5, and Cao 2005). In the spectral calculations, the gravita- 
tional redshift effect is considered, while the relativistic optics 
near the black hole is neglected. This will not affect our final 
results on the XRB, as the inactive galaxies should have ran- 
domly distributed orientations and the stacked spectra would 
not be affected by the relativistic optics. 

The global structure of a RIAF surrounding a black hole 
spinning at a with mass Mbh can be calculated, if some pa- 
rameters, m, a, (3 and 6, are specified. The parameter S 
describes the fraction of the viscously dissipated energy di- 
rectly going into electrons in the accretion flow. The three- 
dimensional MHD simulations suggest t hat the viscosity pa- 
rameter a in the accretion flows is ~ .1 (lArm itage 19981). or 
~ 0.05-0.2 (lHawlev & Balbusll2002l) . The critical accretion 
rate m a -i t ~ 0.01 is suggested by different au thors either fr om 
observations or theoretical arguements (see Naray arj2002l for 
a review, and the references therein). Our numerical calcula- 
tions for the global structure of the flows show m cr i t ~ 0.01 for 
a = 0.2. Thus, we adopt a = 0.2 in this work. The parameter 
/3, defined as p m = B 2 /8n = (l-(3)p tot (p to t = p m +p m ), de- 
scribes the magnetic field strength of the accretion flow. This 
parameter is in fact not a free parameter, which is related to 
the viscosity parameter a as /3 c ± (6a- 3) /(4a- 3), as sug- 
gested by the MHD simulations (|Hawlev. Gammie & B albus 
119961: iNarayan. Mahadevan & Ouataerj|1998l) . For a = 0.2. 
j3 ~ 0.8 is required. It was pointed out that a significant frac- 
tion of the viscously dissipated energy could go into electrons 
by magneti c reconnection, if the magnetic fields in the flow 
are strong dBisnovatyi-Kogan & Lovelacd fl997l l2000h . In 
our calculations, we adopt a conventional value of S = 0.1 
in all the calculations. We plot the X-ray spectra of the 
RIAFs accreting at different rates for Schwarzschild holes 
in Fig. [3] The spectra of the accretion flows surround- 
ing the Kerr black holes with a = 0.9 are ploted in Fig. [4] 
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FIG. 3. — The spectra of the RIAFs surrounding Schwarzschild black holes 
with 1O 8 M0 accreting at different rates: m = 1.0 X 10~ 3 (blue lines), 5.0 X 
lCr 4 (red lines), 3.0 X 10~ 4 (green lines), and 1.0 X 10~ 4 (black lines). The 
viscosity parameter a = 0.2, the fraction of the magnetic pressure 1-/3 = 0.2, 
and the fraction of viscously dissipated energy directly heating the electrons 
8 = 0.1, are adopted in the calculations. The solid lines represent the spectra 
of the RIAFs, while the dotted lines represent their bremsstrahlung spectra 
only. 
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FIG. 4. — The spectra of the RIAFs surrounding Kerr black holes (a = 
0.9) with 10 8 Mq accreting at different rates: m = 5.0 X 10 -4 (blue lines), 
3.0 X 10~ 4 (red lines), 2.0 X 10~ 4 (green lines), and 1.0 X 10~ 4 (black lines). 
The same disk parameters are adopted in the calculations as those in Fig. 
fj] The solid lines represent the spectra of the RIAFs, while the dotted lines 
represent their bremsstrahlung spectra only. 



4. CONTRIBUTION TO THE XRB FROM THE RIAFS IN INACTIVE 
GALAXIES 

The contribution of the RIAFs in all inactive galaxies to the 
cosmological XRB can be calculated by 



fx(E) = 



1 



10 8 M, 



€ ct fe)(l+z)/x[(l+z)£] dV 



4Ttdl 



dz 



dz, (12) 



where lx(E) is the spectrum of a RIAF surrounding a 1O 8 M 
black hole accreting at w^ t , which is calculated as described 
in §3. 1 

There is no doubt that the contribution from active galax- 
ies is important to the XRB, as the BeppoSAX observations 
showed that the power-law X-ray spe ctra of bright AGN s 
extends to severa l hundred keV (e.g., Nicastro et al.l l2000h . 
lUeda e t al. (2003) adopted a template intrinsic AGN spec- 
trum, i.e., a power law+an exponential spectrum with cutoff 
energy E c of 500 keV, in their XRB synthesis model calcu- 
lations. They also included the contribution to the XRB by 
the Compton-thick AGNs with logA^j = 24-25 assuming a 
fixed ratio to those with logA^ = 23-24 without considering 
the redshift dependence of the ratio. The effects of Compton 
scattering leading to significant decrease of the emitted flux 
of the nucleus even in the hard X-ray bands have been prop- 
erly considered in their calculations, while the contribution of 
the sources with logA^ > 25 is neglected by assuming that 
all X-rays are absorbed. In this work, we simply adopt the 
results of their XRB synthesis model calculations for AGNs. 
We consider that the XRB consists not only of the emission 
from type I/II bright AGNs (Compton-thin) described by the 
HXLF, but also of the emission from RIAFs in those inac- 
tive galaxies that are not included in this HXLF. All curves of 
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the XRB synthesis for Compton th in/thick AGNs plo ted in the 
figures of this work are taken from lUeda et al.l (120031) . In Figs. 
[5]and[6] we plot the XRBs contributed by both active galaxies 
and inactive galaxies with different parameters adopted, and 
compare them with the observed XRB. Here, the contribu- 
tion from Compton-thick AGNs has not been included, which 
may be unrealistic, because Compton-thick AGNs have al- 
ready been detected, though their number density is still quite 
uncertain. However, these calculations, at least, provide strict 
constraints (upper limits) on the average mass accretion rates 
of inactive galaxies. We define a parameter £ to describe the 
number ratio of Compton-thick AGNs with logA^H = 24-25 
to those with logA^H = 23-24. In Fig. [7] we plot the results 
of the XRB synthesis, in which the contribution of Compton- 
thick AGNs is considered by assuming £ = 0.5 and 1.0, re- 
spectively. For a given £, the average accretion rate m^?L is 
tuned in order not to ov erproduce the observed XRB in 10-30 
keV dUeda et al.ll2003l) . Thus, the constraints (upper limits) 
on the average accretion rates mflf ct for inactive galaxies from 
the cosmological XRB are plotted as functions of £ in Fig. [8] 
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FIG. 5. — The contribution to the XRB from active galaxies and the RI- 
AFs surrounding Schwarzschild black h oles in inactive ga laxies. The dashed 
line is the observed XRB (taken from Ueda et al. 2003). The dot-dashed 
line repres ents the co ntribution by type I/II AGNs (Compton-thin), which is 
taken from Ueda et al. (2003). The lines with different colors represent the 
XRB synthesis model calculations for different values of the average accre- 
tion rates for inactive galaxies: mf™, = 1.0 X 10~ 3 (blue), 5.0 X 10^ (red), 
and 3.0 X 10~ 4 (green). The solid lines represent the XRB contributed by 
type I/II (Compton-thin) active galaxies decribed by the HXLF and the RI- 
AFs in all inactive galaxies, while the dotted lines are only for the contribu- 
tions from the RIAFs in all inactive galaxies. The black hole mass density for 
inactive galaxies /9tJL act (z) used in the XRB synthesis calcualtions is derived 
for p^ 1 = 4.0 X 1O 5 M Mpc~ 3 . 
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FIG. 6. — The same as Fig. \5\ but for the RIAFs surrounding the Kerr 
black holes with a = 0.9 in inactive galaxies. The lines with different colors 
represent different values of the average accretion rates: wjJSl = 5.0 X 10" 4 
(blue), 2.5 X 10^ (red), and 2.0 X 10" 4 (green). 
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FIG. 7. — The same as Fig. [5] but for the Kerr black holes with a = 0.9 
and the contribution of Compton-thick AGNs is included. The thin dash- 
dotted line repres ents the contribution of Compton-thick AGNs (taken from 
lUeda et aT]|2003D . The yellow line represents the synthesis model for type 
I/II+Compton-thick AGNs assuming £ = 1.0, in which the contrib ution from 
the RIAFs in inactive galaxies is not included (taken from Ueda et al. 2003). 
All lines with other colors correspond to the XRB synthesis model calcula- 
tions by assuming £ = 0.5, i.e., the number of Compton-thick AGNs with 
log/Vn = 24-25 is half of those with logWH = 23 — 24, These lines repre- 
sent different values of the average accretion rates: m?™ t = 5.0 X 10" 4 (blue), 
2.0 X 10~ 4 (red), and 1.5 X 10~ 4 (green). The local black hole mass density 
^bh*' = 4-0 x 10 5 Mq Mpc~ 3 is adopted to derive the black hole mass density 
p'^ ct (z) in inactive galaxies. 
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FIG. 8. — The constraints (upper limits) on the average mass accretion rates 
'"inact as f uncl; i ons °f 'he number ratio £ of the Compton-thick AGNs with 
logyVn = 24 — 25 to those with logNn = 23 — 24. The solid line represents the 
result for pjjjj-' 3 ' = 2.5 X 1O 5 M0 Mpc~ 3 , while the dotted line is for pjjjj-' 3 ' = 
4.0 X 1O 5 M Mpc- 3 . 



5. DISCUSSION 

As described in §2, we derive the black hole mass densi- 
ties for active/inactive galaxies as functions of redshift from 
the HXLF assuming that massive black holes grow through 
accretion. Using different measurements on the local black 
hole mass density, the radiative efficiencies e derived are re- 
quired to be in the range of - 0.07-0.12 (e - 0.11-0.12, 
for p^ Rl = 2.5 x 10 5 M Q Mpc" 3 , and e ~ 0.07-0.08, for 
^tocai _ 4 x 1O 5 M Mpc" 3 , see Table 1 for detailed results). 
If part of the local black hole mass density is accumulated 
through the way other than accretion, the derived radiative ef- 
ficiencies should be even higher than the present values. If 
'"inact ^ 3 x 10~ 4 , less than ~ 5 per cent of local black hole 
mass density was accreted in inactive galaxies. 

The black hole mass density p^ ct (z) for inactive galaxies 
is derived in this work without considering Compton-thick 
AGNs. The total black hole mass density pbh(^) derived in 
this work has been affected little by Compton-thick AGNs at 
low redshifts, because it is derived to match the local black 
hole mass density at z = 0. Our calculations may underes- 
timate the black hole mass density p^(z) in active galaxies. 
Thus, the derived black hole mass density p ™ ct { z) for inactive 
galaxies should be overestimated (see Eq. II lb - Fortunately, 
the derived p l ^ ct {z) has been affected very little at low red- 
shifts, because only a small fraction of massive black holes 
are active at z < 1 (see Fig. |2). The contribution to the XRB 
from inactive galaxies is dominated by that from the sources at 
low redshifts, which means that our main conclusions drawn 
from the XRB synthesis calculations will not be altered, even 
if Compton-thick AGNs are included in black hole mass den- 
sity calculations. 

The derived radiative efficiencies e would be higher than 
the present values listed in Table 1, if Compton-thick AGNs 
are included. If the ratio of Compton-thick to Compton-thin 
AGNs is x, a rough estimate gives the radiative efficiencies 
should be 1 +x times of the present values, because only about 
1/(1 +x) of the local black hole mass density was accumu- 
lated in Compton-thin AGNs and inactive galaxies. Thus, 
the present derived radiative efficie ncies e are only the lower 
limits. If x = 0.6 is adopted (e.g., iTamura et al.l 120061) . the 



radiative efficiencies are about 1.6 times of the present val- 
ues, i.e., e - 0.17 for p^f 1 = 2.5 x 10 5 M o Mpc" 3 , or - 0.1 1 
for p^" 1 = 4.0 x 10 5 M© Mpc" 3 , whic h implies that mos t mas- 
sive black holes are spinning rapidly. lElvis et aTl (120021) sug- 
gested that a high average radiative efficiency e > 0.15 is re- 
quired from the XRB, and they concluded that m ost massive 
black holes must be spinning rapidly. Recently, IWang et all 
(2006) have estimated the average radiative efficiencies of a 
large sample of quasars selected from the Sloan Digital Sky 
Survey (SDSS), by combining their luminosity function and 
their black hole mass function, and found that quasars have 
average radiative efficiencies of 0.3-0.35, which implies that 
the black holes are spinning very rapidly. Our present result s 
seem to be consistent with those derived by lElvis et al.l (|2002l) . 
if suitable number of Compton-thick AGNs is considered. 

iMcLure & Dunlopl d2004l) estimated that the average accre- 
tion rate mf" varies from 0.1 at z ~ 0.2 to 0.4 at z ~ 2 from 
a large sample of SDSS quasars. The average accretion rate 
mil" = 0-1 ls adopted in this work to estimate the inactive 
galaxy black hole mass density p^ cl (z). We find that the 
value of ff^c? 1 ma y a ff ect the resulting p^ ct (z) very slightly 
at low redshifts, because p^ ct (z) Pbhfe) exce Pt at high red- 
shifts (z > 2, see Fig. [2]). The derived black hole mass density 
Pbh Ct (z) for inactive galaxies is mostly governed by the value 
of local black hole mass density p 1 ^^ at low redshifts. The 
XRB is mainly contributed by the sources at low redshifts, 
which implies that the value m^ r only affect our calculations 
of the XRB very little. 

A small fraction of the sources described by the HXLF with 
very high hard X-r ay luminosities ~ io 47 " 48 ergs s -1 are prob- 
ably blazars (e.g.. lUeda et al.ll2003h . For those blazars, the 
beamed X-ray emission from their jets may dominate over 
that from the accretion disks, which may lead to over-estimate 
of the black hole mass density accumulated through accre- 
tion. We perform the calculations of the black hole mass den- 
sities only for the sources with hard X-ray luminosities in the 
range of 10 41 5 - 10 46 ergs s" 1 , and find that the fraction of to- 
tal black hole mass density contributed by these X-ray lumi- 
nous sources with L%[ 10 keV > 10 46 ergs s -1 can be neglected 
(<!%)■ 

We calculate the spectra of the RIAFs surrounding 
Schwarzschild black holes (a = 0), and the Kerr black holes 
with a = 0.9 using reasonable parameters for the RIAFs (see 
discussion in §3), respectively. For the RIAFs surrounding 
Schwarzschild black holes, we find that the X-ray spectra are 
dominated by the bremsstrahlung emission if the accretion 
rates are low (see Fig. [3]), while the Comptonization is im- 
portant for the RIAFs surrounding Kerr black holes (see Fig. 
ID). For a Kerr black hole, the accretion flow extends to a 
smaller radius and the gas density is higher at its inner edge 
compared with that for a Schwarzschild black hole, if they are 
accreting at the same rate. Thus, the Comptonization becomes 
significant for the RIAFs surrounding the Kerr black holes. 
Our calculations show that the peak energy of the spectrum 
decreases with increasing the accretion rate wtjJJ^. For the ac- 
cretion flow with a higher /«f™ t , the density at its inner edge 
is also higher that leads to a lower turnover frequency for the 
synchrotron radiation from the accretion flow, and then leads 
to a lower energy peak produced by the Comptonization, be- 
cause the electron temperatures of the accretion flows at their 
inner edges are roughly similar for different values of MtjJJ^. 

To compare our calculations of the XRB contributed by ac- 
tive/inactive galaxies with the observed XRB, we find that the 
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average accretion rate of inactive galaxies ffij^ should be less 
than 0.001, otherwise, the calculated XRB will surpass the ob- 
served XRB. Our calculations show that m^ t < 2-3 x 10" 4 
is required by the observed XRB (see Figs. [5]and|6]l. It is in- 
teresting to find that the spectral shape of the observed XRB 
(peaked at ~ 30 keV) cannot be well fitted by using any value 
of m^L for Schwarzschild black holes, while both the ob- 
served intensity and the spectral shape of the XRB can be well 
reproduced if rhf™* t — 2.5 x 10~ 4 is adopted for the Kerr black 
holes with a = 0.9 (see Fig. ©for p^ al = 4.0 x 10 5 M o Mpc" 3 ). 
This implies that most massive black holes are spinning 
rapidly, which seems to be consistent with the relatively high 
average radiative efficiencies derived from the XLF (see dis- 
cussion in the second paragraph of this section). A slightly 
higher average accretion rate m a ™ t ~ 3.0 x 10~ 4 is required 
for the case of p[, h cal = 2.5 x 10 5 M Q Mpc" 3 . Note that the 

best-fitted value m^lt c ' oes not depend linearly on 1 / /o}, h cal , 
because the luminosity h as a non-linear de pendence of accre- 
tion rate for RIAFs (e.g. JWu & Caol2006l) . Here, we have ne- 
glected the contribution of the Compton-thick AGNs, as it is 
still unclea r how many Compton-thick AGNs are in the uni- 
verse (e.g., lUeda et al]|2003l) . Thus, these results should be 
only regarded as the upper limits on accretion rates in inactive 
galaxies. We further use a parameter £ to describe the number 
ratio of Compton-thick AGNs with logA^e = 24-25 to those 
with logA^H = 23-24. For the case of £ = 0.5, our XRB synthe- 
sis calculation shows that the XRB contributed by active and 
inactive galaxies (type I/II Compton-thin AGNs +Compton- 
thick AGNs+RIAFs in inactive galaxies) can match the ob- 
served XRB for mass accretion rates rn^L ~ 2.0 x 10 -4 (see 
Fig. Hfor the case of p^ 1 = 4.0 x 10 5 M© Mpc" 3 ). A slightly 
higher accretion rate 2.5 x 10~ 4 is required for the case of 
pioc-A = 2.5 x 10 5 M Q Mpc" 3 . If the contribution from the RI- 
AFs in inactive galaxies is not considered, the same number 
of the Compton-thick AGNs with log.% = 24-25 as those 
with logA^e = 23-24 is required to reproduce the observed 
XRB (see the yellow line in Fig. [7). In Fig. [8] we find 
that the contribution to the XRB from inactive galaxies can 
be neglected compared with that from Compton-thick AGNs, 
if m^lt ^ 10~ 4 . Corresponding to such constraints on the av- 
erage accretion rates for inactive galaxies, we find that less 
than ^5 per cent, or even less than ~2 per cent if Compton- 
thick AGNs are considered, of the local massive black hole 
mass density was accreted during the radiatively inefficient 
accretion phases (see Table 1). The constraints of the XRB 
can only give upper limits on /«?™ t f° r inactive galaxies, and 
we cannot rule out the possibility that inactive galaxies are 
accreting at rates significantly lower than 10~ 4 . This may be 
resolved by deep surveys on faint X-ray sources in hard X-ray 
bands, which is beyond the scope of this paper. 

The detaile d AD AF spectral calculations by 
iDi Matteo et al.l ( 1 19991) showed that many sources at 
redshift z ~ 2 - 3 with ADAFs accreting at the rates close 
to the critical value would be required to reproduce the 
observed spectral sha pe of the X RB with an energy peak at 
~ 30 keV. However, ICaol (l2005l) 's spectral calculations for 
RIAFs accreting at the critical value showed higher energy 
peaks at several hundred keV, which are unable to reproduce 
the observed spectral shape of th e XRB (see Fig. 2 in Cao, 
2005). The reason is that ICaol d2005l) ado pted a larger S 
than that used by Di Matteo et al. (1999) for traditional 
ADAFs, which lead s to higher electron tem peratures of the 
accretion flows (see Cao 2005, for details). IDi Matteo et al.l 



( 1999)'s calculations are limited to the spectral shape, which 
have not included the contribution of the ADAFs in all 
galaxies to the XRB quantitatively, and the contribution of 
normal bright AGNs has not been considered either. Our 
present calculations show that the spectra of the RIAFs 
surrounding Kerr black holes have peaks at around 10-100 
keV, which depend sensitively on accretion rates (see Fig. 
4). We find that our model calculations can fit both the 
intensity and the spectral shape of the XRB very well, if 
the RIAFs surrounding the Kerr black holes with a ~ 0.9 
in inactive galaxies are accreting at mf^ t ~ 1 - 3 x 10 -4 . 
Cao (2005) calculated the contribution to the XRB from 
the RIAFs accreting at the critical rate, and found that the 
timescale for the sources accreting at the critical rate should 
be much shorter than the bright AGN lifetime. In that work, 
it was found that the spectral shape of the XRB cannot be 
reproduced even the RIAFs' contribution is included. Our 
present calculations imply the contribution from the RIAFs 
in the inactive galaxies accreting at ~ 1 — 3 X 10~ 4 should 
dominate over that from those accreting at the critical value. 
It means that the timescale of AGNs accreting at aroun d the 
critical value should be even shorter than that given bvlCaol 
J20051) . 

In this work, we have to assume the black holes in all in- 
active galaxies are accreting at an average rate »i a ™ t , because 
we are almost ignorant of how the accretion rate evolves with 
time in active/inactive galaxies. A more realistic possibility 
may be that the massive black holes in inactive galaxies are 
accreting at different rates from far below ~ 10~ 5 up to the 
critical value of ~ 0.01, of which the stacked spectra (com- 
bined with the contribution from the type I/II AGNs) may re- 
produce the observed XRB. The spectrum of the RIAF ac- 
creting at the critical value is peaked at ~ lOO keV, d ue to the 
saturated Comptonization (see Fig. 1 in Cao 2005). These 
RIAFs are X-ray luminous, and therefore much more inactive 
galaxies accreting at very low rates (<C rn^Zt) are required 
to reproduce the intensity of the XRB. However, the spectral 
shape of the XRB with an energy peak ~ 30 keV cannot be 
fitted in this way, because the energy peaks of the RIAF spec- 
tra are around several hundred keV while they are accreting 
at m < 10~ 4 (see Figs. [3] and |4]i . It means the contribution of 
the RIAFs accreting at around the critical value to the XRB 
should be unimportant compared with other inactive galaxies, 
which, again, suggests that the accretion rates of the RIAFs 
in inactive galaxies should decreases very rapidly after the ac- 
cretion mode transition. It provides evidence of the feedback 
on the circumnuclear gases from the nuclear radiation, which 
is qualitati vely consistent with the numerical simulations (see 
Fig. 2 in lDi Matteo et al.ll2005l) . 

The RIAF may have winds, and a power-law radius- 
dependent accretion rate, rh(r) oc r~ Pw (p w > 0), is 
assumed, though the detail ed physics is still unclear 
(Blandford& Begelman 1999). As most X-ray emission is 
from the inner region of the accretion flow near the black hole, 
the X-ray spectrum of a RIAF is similar to a RIAF with winds, 
provided they have the same accretion rate at their inner edge. 
In this work, we are investigating the growth of black holes, 
so we focus on the accretion rate at the inner edge of the ac- 
cretion flow. Thus, the conclusions will not be altered, even if 
winds are present in the RIAFs. 



This work is supported by the National Science Fund for 
Distinguished Young Scholars (grant 10325314), and the 



Cao 



TABLE 1 

The summary of different models 



Model 


"'inact 


€ 


^locaia 


Pb B h (0)/p b ° h cal 


e 


Jocal a 
Pbh 


Pbb(0)/P b °b al 


A 


1.0 x l(r 3 


0.122 


2.5 x 10 5 


0.162 


0.080 


4.0 X 10 5 


0.162 


B 


3.0 x ur 4 


0.110 


2.5 x 10 5 


0.051 


0.071 


4.0 X 10 5 


0.051 


C 


2.0 x ur 4 


0.108 


2.5 x 10 5 


0.034 


0.070 


4.0 x 10 5 


0.035 


D 


1.0 x ur 4 


0.106 


2.5 x 10 5 


0.017 


0.069 


4.0 x 10 5 


0.017 


E 


0.0 


0.105 


2.5 x 10 5 


0.0 


0.068 


4.0 x 10 5 


0.0 



a in units of Mq Mpc 3 . 



NSFC (grant 10333020). 
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